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A METHOD, AM INSTALLATION, AND A COMPUTER PROGRAM FOR 
ESTIMATING THE INITJJ\L SIZE OF A POPULATION OF NUCLEIC 
ACIDS, IN PARTICULAR BY PCR 
The present invention relates to estimating the 
5 ■ initial size of a population of interest in a sample 
subjected to a succession of amplification reactions. 
BACKGROUND OF THE INVENTION 
The present invention finds a particularly 
advantageous, but non-limiting, application in 
10 determining an initial quantity of nucleic acids in a 
sample subjected to a polymerase chain reaction (PCR) in 
real time. A technique of this type, known as "PCR 
quantification" , is used in particular for evaluating the 
number of copies of pathogenic agents (e.g. of the human 
15 immunodeficiency virus (HIV) ) in a sample of body fluids 
taken from a patient, typically in the context of a' 
medical checkup. 

Reference is made to Figure 1 for a brief 
description of the diagrammatic appearance of a real time 
20 PCR amplification curve with PCR cycle index, numbers 
plotted along the abscissa and, in the example shown, 
with quantities of . fluorescence emitted (in arbitrary 
units) as measured for each PCR cycle plotted up the 
ordinate. For each PCR cycle, it should b,e understood 
25 that the sample is subjected to temperature variations 
enabling DNA polymerase to amplify nucleic acids and 
enabling the corresponding PCR products to be detected by 
fluorescent molecules. By plotting the measured 

fluorescence F,^ as a function of PCR cycle number n, 
30 variation is obtained of the type shown in Figure 1, and 
comprises at least: 

• a first portion BN where fluorescence measurements 
coincide substantially with the background noise of the 
apparatus for measuring fluorescence; 
35 • a second portion EXP in which the measured 

quantities of fluorescence increase in substantially 
exponential manner; 
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• a third portion LIN in which the increase in the 
measured quantities of fluorescence is significantly 
attenuated and behaves, overall in substantially linear 
manner; and 

5 ■ a fourth portion PLA in which fluorescence 

measurements reach a plateau stage. 

It should be observed that for the initial PGR 
cycles (first and second portions), the population of 
interest increases in substantially exponential manner, 

10 whereas for the following cycles (third and fourth 
portions) , other phenomena come into competition with 
growth in the population of interest, so that said growth 
is then damped up to the plateau stage PLA. 

The document "Mathematics of quantitative kinetic 

15 PGR and the application of standard curves" by 
R.G. Rutledge and C. Gote, published in Nucleic Acids- 
Research, 2003, Vol. 31, No. 16, discloses a method of 
estimating the unknown initial quantity of nucleic acids 
in a sample of interest by means of PGR. That method 

20 consists in using a plurality of samples havijig known 
initial quantities of nucleic acids, referred to as 
"standards", in order to determine by interpolation the 
initial quantity of nucleic acids present in the sample 
of interest. 

25 In general, the greater the initial quantity of 

nucleic acids in a sample, the sooner a detectable 
quantity of PGR product is obtained, i.e. the sooner a 
detectable quantity of emitted fluorescence is obtained. 
With reference to Figure 2, relating to the prior art, it 

30 will be understood that the initial population in the 
standard St 1 is _ greater th^^ ,.the .standard St2.. 

which is greater than that in the standard St3, etc., 
since the cycle Ctl for the standard Stl occurs before 
the corresponding cycle Ct2 for the standard St2, which 

35 occurs before the cycle Ct3 for the standard St3, etc. 

Thus, such a CT cycle, corresponding to the cycle at 
which the fluorescence measurements reach a fluorescence 
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threshold " THR (a^ shown in Figure 2) , sets at an 
arbitrary level (typically below the background noise) , 
and acts as a parameter representative of the initial 
size Nq of a population of. nucleic acids subjected to the 
5 PGR cycles. Use has been made of this observation in the 
above-cited prior art to establish a relationship of the 
kind shown in Figure 3 between cycle numbers Ctl, Ct2, 
CtS, ' Ct4 for a plurality of standards having known 
initial populations, and their initial populations Nq^, 

10 Nq^, Nq^, No'*. Thus, by plotting the cycles Ctl, Ct2, Ct3, 
Ct4, etc. up the ordinates and the logarithm of the 
initial population sizes Nq^, No^f Nq^, Nq* along the 
abscissa, a regression slope REG is obtained. On this 
regression slope PEG, the cycle Ctint detected for the 

15 sample of interest is plotted (dashed-line arrow Fl) . By 
interpolation on the regression slope REG (dashed-line 
arrow F2), the initial population size Nq^"*^ is then 
determined for the sample of interest. 

Although that method is in widespread use, it 

20 nevertheless presents some drawbacks. 

Firstly, it requires the use of a plurality of 
standard samples having respective known initial 
populations. 

Secondly, the method depends on the judgment of the 
25 user, since the fluorescence threshold value, as selected 
by the user, has a direct influence on the values of the 
Ct cycles in the amplification curves, and consequently 
on. the estimated values for the initial population size 
in the sample of interest. The threshold value also has 
30 an impact on the accuracy of the result, since accuracy 
is generally better if the threshold is selected to lie 
in the exponential growth stage EXP of the amplification 
curve. Nevertheless, in practice, it is difficult for 
the user- to know whether the fluorescence threshold level 
35 THR that .has been set does indeed correspond to the 
exponential stage of the curves, and does so for all of 
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the samples- (the standard samples and the sample of 
interest) . 

Finally, the method assumes without any verif ica tiLon 
that the population has the same amplification yield in 
5 the sample of interest and in all of the standard 
samples. Thus, if the sample of interest contains PGR 
inhibitors, as is typically the case, then its result 
will be falsely lowered. 

It should thus be understood that the prior ^rt 

10 technique depends on the fluorescence threshold THR as 
defined by the user. The value selected has an influence 
on the values of the Ct cycles and consequently on 
determining the initial quantity in the sample of 
interest. That is one of the reasons why a large amouint 

15 of work has recently been undertaken to automate Ct cyc^le 
detection and make it reliable. 
OBJECTS AND SUMMARY OF THE INVENTION 

The present invention seeks to improve the situati^on 
by proposing an approach that is completely different. 

20 Firstly, the invention provides a method, the methxod 

being implemented by computer means for quantifying in 
absolute and/or relative manner an initial population of 
nucleic acids in a sample of interest. . The sample is 
subjected to a succession of applications of a reacti_on 

25 for amplifying the population of interest. In very 
general manner, this amplification may be undertaken by 
implementing successive PGR cycles, however any otht.er 
amplification technique could also be used. Above al.1, 
it should be understood that the amplification needs 

30 merely to be defined by a reaction yield, as described 
below. During these successive amplification operatiorLS, 
experimental measurements are taken that are 
representative of a current population size, at least in 
the sample of interest . It will be understood that one 

35 or more measurements can be taken after or during each 
amplification reaction without loss of generality. 
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In a presently preferred definition of the 
invention, the method in the meaning of the invention 
comprises the following steps: 

a) providing a model of the yield of the 
5 amplification reaction as a function of the succession odf 

amplifications, said model comprising: 

• a substantially constant stage for a first 
portion of the applications of the amplification 
reaction; and 

10 • a non-constant stage for a second portion of 

the applications of the amplification reaction; 

the first and second portions being united by 3 
changeover region in which yield changes over between the 
constant and non-constant stages,, said region having an 

15 amplification index corresponding substantially to the 
changeover; 

b) using the yield model to express a relationship, 
involving at least the changeover index and a parameter- 
representative of the initial population size in the 

20 sample of interest; 

c) determining at least the changeover index hy 
comparison with the experimental measurements; and, in a. 
subsequent or immediately following step d) deducing- 
therefrom the initial population size in the sample of 

25 interest, 

BRIEF DESCRIPTION OF THE DRAWINGS 
Other advantages and characteristics of the 
invention appear on reading the following detailed, 
description of an implementation given below by way of 
30 example with reference to the accompanying figures, ia 
which: 

• Figure 1 relates to the prior art and represents 
variation in the measured - quantity of fluorescence as a 
function of the number of PGR cycles, as described above; 
35 • Figure 2 relates to the prior art and is 

representative of the increasing quantities of 
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fluorescence that are emitted as a function of the number 
of PGR cycles, as described above; 

Figure 3 represents an interpolation method for 
determining the initial quantity of the population of 
5 interest in the sample of interest using a method known 
in the prior artr and described above; 

■ Figure 4A is a diagram showing variation in the 
above-described experimental measurements as a function 
of the succession of amplifications applied to the sample 
10 of interest; 

• Figure 4B is a diagram showing variation in the 
yield of the amplification reaction, obtained from 
experimental measurements, as a function of the 
succession of amplifications applied to the sample of 

15 interest; 

• Figure 5 plots a regression relationship between 
the yield changeover indices that occur at the changeover 
between the constant stage and the non- constant stage, 
and the logarithms of the initial populations for 

20 standard samples and for the sample of interest,, for use 
in a first implementation; 

• Figure 6A shows typical variation in the measured 
quantity of ' fluorescence after it has been adjusted by 
taking account of background noise specific to the 

25 measurements, and plotted as a function of the number n 
of PGR reaction cycles; 

• Figure 6B shows the variation in the effectiveness 
of the PGR shown in Figure 6A as a function of the number 
n of cycles; 

30 • Figure 7 compares the experimental variation in 

emitted fluorescence as s^pwn in Figure _6A ...with _ the. 

results obtained by applying the emitted fluorescence 
model obtained by including an effectiveness model in a 
second implementation; 

35 • Figure 8 is .a comparison between the variation in 

the effectiveness of Figure 6B and the application of an 
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effectiveness model deduced frora the emitted fluorescence 
model of Figure 7; 

• Figure 9 is a flow chart outlining the main steps 
in the method in a proposed implementation of the present 
5 invention; and 

Figure 10 is a diagram of an installation for 
quantifying the initial population of a sample of 
interest . 

MORE DETAILED DESCRIPTION 

10 Reference is made to Figures 4A and 4B for briefly 

describing a few principles of the invention illustrating 
the characteristics of the above method. 

Firstly, it should be understood that Figure 4A 
plots a succession of" experimental measurements F,^ 

15 representative of the current size of a population of 
interest which is being subjected progressively to a 
succession of . amplification reactions, each reaction 
being indexed by an index number n. In the non-limiting 
example described herein, this succession, of reactions 

20 corresponds to a succession of PGR cycles. In this non- 
limiting example, the experimental measurements F^ 
correspond to measured quantities of fluorescence on each 
PGR cycle. Thus, in a quantification method which 
combines the PGR reaction and the fluorescence emitted by 

25 the sample of interest, fluorescent reagents are 
introduced into the sample so that the fluorescence that 
is emitted during a PGR cycle is proportional to the size 
of the nucleic acid population in the sample. Indeed, it 
can be preferable to perform a plurality of measurements 

30 or no measurements at all for certain PGR cycles. 
Furthermore , more generally_,^^ _ the. _mea_sur^^ , method may. 
make use of techniques other than fluorescence, even if 
fluorescence is the method that is pften used for 
quantifying by PGR. Finally, it should be understood 

35 that other amplification techniques could be implemented 
in the context of the present invention, providing it is 
possible to track variation in the yield of the reaction 
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corresponding to the amplification. Since the example 
described below relates preferentially to PGR cycles, 
j^eference is made to the "effectiveness of the PGR" 
written E^^ for each PGR cycle of index Uf in order to 
5 refer to the yield of the amplification reaction. 

As mentioned above with reference to Figure 1, 
Figure 4A mainly comprises two regions in which: 

• during the initial PGR cycles (portion EXP) , the 
population increases substantially exponentially; whereas 

10 ' during the following cycles (the LIN and PLA 

portions), other phenomena come into competition with 
growth of the population of interest, so the growth 
becomes damped. 

The following two assumptions are made: 

15 • the yield of the reaction E^^ is relatively constant 

during the initial cycles over the portion EXP; and 

• after some number of cycles have been performed, 
the yield of the reaction starts to decrease over the 
portions LIN and PLA. 

20 . This decrease in yield may have a variety of 

explanations, in particular a degradation and/or a lack 
of PGR reagents (DNA polymerase, dNTPs, primers, etc.) 
and/or inhibition by the products that are made 
themselves . 

25 "It. is assumed herein that the yield is initially 

constant and that " it subsequently decreases. 
Nevertheless, it should be understood that the invention 
applies more generally to the context of yield: 

• that is initially constant, which corresponds to a 
30 normal situation for growth by amplification; and 

• that is subsequently _ not .^cpnst^^ .or 
increases) which corresponds to a situation that is 
substantially abnormal. 

In the context of reactions for amplifying the 
35 quantity . of nucleic acids, it has been found that the 
yield often changes over from a constant stage to a non- 
constant stage. In the meaning of the invention. 
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advantage ' is taken of this observation to deduce 
therefrom the initial quantity of nucleic acids, as 
described below in detail. Initially, it is merely 
stated that the yield can also change over from a non- 
5 constant stage during early cycles to a subsequent 
constant stage. The present invention is equally 
applicable to such a circumstance. In general, it should 
therefore be understood that in the meaning of the 
invention, a changeover of yield between a constant stage 

10 and a non-constant stage is detected. 

The objective is to find the initial size of the 
population that has been subjected to amplification. 
With reference to Figure 4, it will be understood that 
the measurement Fq representative of this initial 

15 population size, which coincides in practice with the 
measurement background noise BN, cannot be used on its 
own for determining directly the initial population size. 
In the prior art, attempts have been made to quantify 
this initial population size by making use of the 

20 exponential stage, i.e. a stage that occurs typically on 
exiting background noise. A threshold cycle Ct is then 
determined (corresponding to point PA for "prior art") in 
Figure 4A. As mentioned above, in this region 

measurements are often affected by noise and it is 

25 , difficult to determine accurately a threshold cycle Ct 
representative of exiting background noise. 

In a completely different approach, the present 
invention instead makes - use of nearly all of the points 
of the amplification curve in order to determine 

30 accurately a region CHO where the yield changes over 
between, .a _ constant .„ stage_. ..and .a... non-conatant . stage, 
typically in present circumstances between the 
exponential stage EXP and the linear stage LIN, It will 
be understood that measurements are logically less 

35 affected by noise in this region CHO than in the 
background noise exit region since the region CHO occurs 
during later cycles. Furthermore, particularly because 
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of the mathematical properties associated with yield, it 
is shown below that, most advantageously, the number of 
standards that need to be used for quantifying the 
initial size of the population of interest is smaller 
5 than the number of standards used in prior art 
quantification. 

The relationship for associating the changeover 
region CHO with the initial size of the population of 
interest is briefly described below. The yield of an 
10 amplification reaction is given by: 

in which: 

• N„ is the size of the population of interest after 
an amplification of index n in a succession of 

15 amplifications; 

Nj^+;l is the size of the population of interest 
after a following amplification, of index n+1, in the 
above-mentioned succession of amplifications; and 

• En is .the yield of the amplification reaction of 
20 index n in the above-mentioned succession of 

amplifications. 

Reformulating this relationship as a recurrence 
relationship, we obtain: 

Nn+l = (1 + En) (1 + En-i) (1 + En-2) ... d + Eo)No 

25 where Nq is the initial size of the population of 
interest. So long as the yield E^ is constant, it will be 
understood that the above relationship can be written 
more simply as follows: 
Nn+i = No X (1 + Eo)^^i 

30 where the index n+1 has not yet reached the changeover 
region CHO. While the yield is. cpnst_ant ...during . the. 
initial cycles, the following applies: 

En = En-l = En-2 = ... = Eq 

where Eq is the value of the yield during the constant 
35 stage. Nevertheless, when the index n+1 moves into the 
changeover region CHO, the relationship becomes: 
Nn+i = No X (1 + Eo)^EEP X function (Ceep/ n+1) 
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where: 

■ (^EEP " ■'■^ last index of the amplification 

reaction during which the yield is still constant (it 
will thus be understood that the index Cg^p its elf 
5 represents the changeover index proper between the 
exponential stage and the linear stage); and 

the term function {C^gp^ n+1) is a particular 
function characterizing the non-constant stage of the 
yield and that depends at least on the changeover index 
10 Ceep and on the current amplification index n+1. 

It- can thus be seen how it is possible to associate 
the changeover index C^Ep and the initial size Nq of the 
population of interest. At this stage it can be 
understood that steps a) and b) of the above-defined 
15 method have already been implemented, 

A first implementation consists in determining the 
changeover index Cg^p experimentally and' in correlating it . 
with the initial size by regression by using a plurality 
of standard samples that are subjected to the same 
20 amplification treatment as the sample of interest. Under 
such circumstances^ it will be understood that steps b) 
and c) of the above-defined method are merely 
interchanged since initially the changeover index Ceep 
(step c) ) is determined experimentally, and subsequently 
25 the relationship between the index C^^p and the initial 
size Nq (step b) ) is determined in order to end up with 
the initial size Nq (step d) ) . 

Before describing all of these steps in detail in 
the meaning of the first implementation, a method is 
30 described for determining the index Cggp on the basis of 

.experimental.. .measuxements partl-cular.,----it will be- 

understood that this method of determining the index C^gp 
experimentally can be applied to another implementation 
that is different from the above-mentioned first 
35 implementation. 

Returning to the relationship between the 
effectiveness E„ of a given cycle n and the current size 
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of the population of interest in the same cycle and in 
a subsequent cycle N^^+i, the effectiveness of the 
amplification can be expressed as follows: 

5 In certain circumstances^ in particular when there 

is no need to take account of background noise BN in the 
measurements, it is possible to a first approximation to 
assume that the measurements are substantially 
proportional to the current size of the population of 

10 interest. Nevertheless, in practice, account will more 
often be taken of measurement drift, -with corrected 
experimental measurements F\ being determined on the 
basis of direct measurements as shown in Figure 4h. 

A prior step of processing the experimental 

15 measurements is preferably applied, this step 

consisting in subtracting the background noise BN and 
subsequently in introducing compensation to take account 
of a non-zero measurement s representative of the initial 
population size. In the example shown in Figure 4A, the 

2 0 variation in the background noise BN as a function of the 
index n can be represented by a linear function since 
tests have shown that a linear model is satisfactory for 
fluorescence measurements in PGR- Nevertheless, in 
certain circumstances it may be preferable to use an 

25 exponentially-varying model. In any event, a model is 
applied that complies best with variation in the 
background noise BN as given typically by the initial 
measurement points. Thereafter, the selected model for 
variation in background noise BN is subtracted from all 
' 30 of the experimental measurement values F^. By applying 
this step, it will .be understood,_that. „the. theoretical,, 
fluorescence measurement Fq is reduced to a measurement 
value of zero, corresponding to an initial population 
size No of zero, which is not representative of physical 

'35 reality. Consequently, it is advantageous to apply 
compensation for this correction as follows: 
■ F'„ = F„ - BN + s 
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where: 

■ the term corresponds to a corrected measurement 
for a current index n; 

• the term corresponds to the raw experimental 
5 measurement at said current index n; 

the term BN corresponds to the value for the 
background noise as modeled for the current index n; and 

• s is the corresponding compensation term which is 
assumed to be constant in the example being described and 

10 which directly represents the initial population size N^. 

Although these steps of correcting for background 
■ noise are very advantageous in determining the changeover 
index C^^^, they may also be applied to any determination 
and quantification of the initial population size Nq 

15 whenever background noise is likely to falsify 
measurement of said population size Nq. In this respect, 
these steps may constitute the subject matter of separate 
protection, where appropriate. 

The corrected measurements F'^ as obtained in this 

20 way are advantageously proportional to the . current 

population sizes in the samples of interest, such that 

the yield can now be expressed directly as a function 

of measurement values (corrected as described above), by 

the following relationship: 

25 = (F^,a/F'n) - 1 

Thus, from the experimental measurements F^ of 

Figure 4A, corrected measurements F'^ ^re obtained from 

which there is subsequently determined the variation in 

the effectiveness ^ function of the succession of 

30 indices n, as shown in Figure 4B. 

In short, the experimental __it^^^ ,.are 
expressed in the form of an experimental variation in the 
' effectiveness E„ of the kind shown in Figure 4B as a 
function of the succession of amplifications n. This 

35 provides an experimentally-determined' variation for the 
yield comprising: 
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a" perceptibly noisy first region for low 
amplification indices n (specifically prior to the cycle 
CG in the example of Figure 4B) ; and 

• followed by a second region exhibiting less noise 
5 for higher amplification indices (at least after the 
changeover region CHO) . 

At least in the most usual circumstance of 
amplification by PGR and measurement by fluorescence, the 
non-constant stage of yield is decreasing and corresponds 
10 to said second region presenting little noise (as shown 
in Figure 4B) . Specifically for the purpose of 

eliminating measurement points that run the risk of 
falsifying results when selecting a model to apply to the 
variation in yield: 
15 . a crude value Eq is estimated for the constant 

yield stage; and 

particularly when searching for the changeover 
index C^ep' least some of the measurements in the less 
noisy second region are ignored for which the estimated 

20 yield is less than a threshold value, e.g. corresponding 
to some fraction of the constant stage Eq. 

These points NEG (Figure 4) that are eliminated are 
typically those that correspond to very high 
amplification indices n and that might no longer satisfy 

25 the model for effectiveness which is selected 
substantially around the changeover region CHO. By way 
of example, in order to eliminate them, an average is 
evaluated for the constant stage of yield Eq, typically 
for the initial indices n. Thereafter, a threshold value 

30 is selected that corresponds to a fraction of the average 
_ found for_the_ constant stage Eq/. e.g.. 10%.._ . Thereafter, 
starting from the highest indices n, all measurement 
points NEG. of measured yield lower than or equal to said 
threshold value are eliminated. This step, which is most 

35 advantageous for detecting the index C^^^, can 
nevertheless be applied to any determination based on 
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yield E„/ and may constitute the subject matter of 
separate protection, where appropriate. 

When yield presents a non-constant stage in whi ch 
yield is decreasing and which follows a constant stage, 
5 as shown in Figure 4B, the changeover region CHO is 
identified by working in the direction of decreasing 
index numbers n, starting, from the less noisy second 
region, and by detecting a coarse index CG for which the 
yield passes a predetermined value. Thus, with reference 
10 to Figure 4B, going in the direction of decreasing index 
numbers n so as to rise towards the changeover region 
CHO, the yield associated with each index number is 
evaluated. For the first measurement point of yield that 
is significantly greater than the above-mentioned 
15 predetermined value, it is considered that the above- 
mentioned coarse index CG has been detected and 
corresponds to the index of the measurement point. 

As described below in a subsequent implementation, 
it is possible for each measurement point to model the 
20 variation in its yield as though said set point itself 
corresponds to the changeover index C^^^. In that 
implementation, if the constant yield stage Eq is 
estimated, and if the estimated value exceeds the above- 
mentioned predetermined value, then the point is 
25 considered as corresponding to the coarse index CG. 

In general, a maximum yield has a value of 1 so it 
is possible to select the above-mentioned predetermined 
value as being equal to 1. Nevertheless, this can loe 
varied, and, for example, provision can be made to set 
30 the predetermined value as corresponding to the mean 
yield Eq as evaluated over the initial reaction cycles - 

Thereafter, the estimate of the value of the 
amplification index . C^^p in the changeover region is 
refined, which value may advantageously be a fraction, by 
35 working in the .direction of increasing amplification 
index numbers, starting from the coarse index CG, and by 
detecting, an amplification index for' which the yield is 



wo 2006/048337 



PCT/EP2005/013021 



16 

approximately equal to the above-mentioned predetermined 
value. Thus, referring again to Figure 4B, in order to 
refine the search for the changeover index C^ep after the 
coarse index CG has been determined, a search is made 
5 downwards starting from the coarse index CG and going in 
the direction of increasing index number n, in steps of a 
size smaller than one whole index, and the abscissa value 
is determined, e.g. by interpolation, at which the 
predetermined value is crossed. Typically, so long as 

10 the constant stage value remains greater than 1, the 
search is continued in the direction of increasing number' 
n, and the index {C^^^ - I) preceding the changeover is 
determined as soon as the constant value Eg is equal to or 
very close to 1. That is why it is appropriate to select 

15 a search step size corresponding to a fraction of the 
index, for example 10% of one cycle n. 

In the above-mentioned first implementation, a 
plurality of standard samples are provided having 
respective known initial population sizes, and the 

20 succession of amplifications is applied thereto under 
substantially the same conditions as for the sample of 
interest. Their respective changeover indices are 
determined in accordance with above-described steps a), 
b) , and c) . In step d) : 

25 • a dependency relationship is established between 

the, initial population sizes of the standard samples Nq^^ 
and their indices 0^^?^^; and 

• after determining the index C^gp for the sample of 
interest, the initial size of the . population of interest 

30 Nq is determined ' by interpolation on that dependency 
relationship . . 

Thus, with reference to Figure 5, a dependency 
relationship can be established between the changeover 
cycles Ceep^ ^be?^ standards and their initial 

35 concentrations Nq^, Nq^ (actually the logarithms thereof), 
e.g. by regression. By measuring the changeover index 
CnjEp^"*^ for the sample of interest and by plotting its 
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value on the regression slope of Figure 5, the inditial 
concentration Nq-^"*^ in the sample of interest is obtained 
by interpolation. 

This first implementation is thus quite similar to 
5 that of the prior art described with reference to 
Figure 3. Nevertheless, it should not be forgotten that 
the changeover index Cggp on which this first 
implementation relies does not correspond in any way to 
the threshold cycle Ct of the prior art, 
10 In an approach that is significantly different from 

this first implementation: 

• in step b) , use is made of the yield modeJ to 
express variation that is parameterized as a function of 
the succession of amplifications, said variation raa-king 
15 use of at least one parameter representing the changeover 
index Ceep/ and 

in step c) , at least said parameter represen. ting 
the changeover index C^^p is determined by comparison with 
the experimental measurements. 
20 In a second implementation, this parameterized 

variation is representative of the current population 
size in the sample of interest. 

Typically, this parameterized variation can be drawn 
from an expression of the type given above: 
25 Nn+i = No X (1 + Eo)^EEP X function ( Cbep^ n+1) 

Thus, in addition to a parameter representing the 
changeover index C^Ep, this variation makes use of a 
parameter representative of the initial population sizie Nq 
in the sample of interest. 
30 Thereafter, in steps c) and d) of this seoond 

.implementation,,. ..these.. ...two ..parameters.. .Ceep- • ^and . Nq. -axe- 
determined substantially together. 

Previously, in step a), it is necessary to determine 
a model for the above-mentioned function function (CgEp^ 
35 n+1) . 

Usually, for PGR quantification, a model is selected 
for the non-constant stage of the yield corresponding to 
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a decreasing exponential having a decrease parameter P 
which is described in greater detail below. This 
decrease parameter P is then determined in step c) , at 
least with the changeover index C^^^, by comparison with 
5 the experimental measurements. 

Thus, in this second implementation, once the yield 
model En has been selected, it is applied to the general 
expression for. the current population size given by the 
above relationship. This provides a model for variation 

10 in the current population size N^- 

Nevertheless, unless the experimental measurements 
give the value for the current population size directly 
(which is rarely true in practice at present), it is 
appropriate subsequently to model the experimental 

15 measurements themselves, taking account of the 

subtracted background noise and the subsequent 
compensation 8 as described above. 

Thus, in a presently preferred implementation, the 
above-mentioned parameterized variation: 

20 • is representative of experimental measurements; 

and 

includes a parameter corresponding to a 
measurement value Fq representative of the initial 
population size. 
25 Thereafter, the measured value of the initial 

population size Fq is determined by comparing said 
parameterized variation F^ with the .experimental 
measurements . 

In order to perform this comparison, it is possible, 
30 for example, to adjust the parameters Fq, Eq, Ceep and the 
decrease parameter p _ in the model of jyie measureme^^ F^ 
by using statistical correlations (typically the least 
squares method) applied to the raw experimental 
measurements. An example implementation is described in 

35 detail below. 

Initially, variation is obtained for the measured 
. and adjusted quantity of fluorescence as a function of 
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the number of PGR cycles that have been applied, as shown 
for example in Figure 6A. This figure shows the 
■ amplification curve for a sample of interest containing 
nucleic acids, in this case a fragment of DNA having an 
5 initial quantity -of 100,000 copies, marked by the 
SYBRGREEN intercalant during the PGR reaction which is 
performed on the I-Cycler IQ® apparatus from the supplier 
BI-RAD®. 

In the example described, it will be understood that 
10 the amplification reaction is a PGR reaction in real 
time. The experimental measurement represents quantities 
of emitted fluorescence. 

The fluorescence of cycle n after adjustment for 
background noise, as described above, is written below. 
15 The theoretical initial fluorescence before the first 
cycle is written Fq. The effectiveness of the PGR in 
cycle n is written E^,. The total number of cycles 
performed during the PGR reaction is written N. 

By assumption, the fluorescence measured on each 
20 cycle n of the PGR reaction cycle is defined by: - 

^n^i ^ Fnd + Ej for all n e {o, 1, 2,...,N -l} (1) 
with 0 ^ E„ ^ 1. 

The effectiveness of the reaction on each cycle n is 
calculated as follows: 

25 E - 1 « — - 1 for all n e (0, 1, 2,..,,N -l} (2) 

" N F 

It should be observed that equation (1) is assumed 
to be true for n=0. Nevertheless, by definition, the 
initial fluorescence Fq is unknown. It is therefore not 
possible to calculate the effectiveness on the first 
-30- ■ cycle Eq directly f ox* formula- (2) ; 

Figure 6B shows the effectiveness of the PGR 
reaction as approximated by formula (2) and on the basis 
of the adjusted variation in fluorescence of Figure 6A, 
as a function of cycle number n. 
35 The following assumptions are preferably made: 
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• the effectiveness of the reaction is relatively 
constant during the initial cycles; and 

after a certain number of cycles have been 
performed, the effectiveness of the reaction decreases. 
5 Figure 6B confirms the second assumption since it 

can be seen that effectiveness decreases as from cycle 
n=17. However, the measured effectivenesses in cycles 1 
to 16 are very noisy, which makes it difficult to verify 
the first assumption graphically. 
10 Nevertheless, it is preferable to assume that 

variation in effectiveness obeys a model of the type 
including:. 

• a first stage that is constant between the first 
PGR cycle and the cycle (C^^^-l) preceding the changeover 

15 cycle written Ceep; and 

• a second stage in which it decreases for cycles of 
numbers greater than or equal to cycle (CgEp-l) • 

The cycle (Ceep~1) thus represents the last cycle 
(which may be a fraction) for which effectiveness 
20 continues to be constant. 

It is then proposed to model the effectiveness .of 
the reaction as follows: 

for 0 < n < (Ceep -1) 
(l-fEor^^'*'^"-^"'^'* - 1 for (Ceep-1) < n < (N-1) 

where Eq and P are real parameters which are estimated 
25 using the amplification curve of Figure 6A, or using the 
effectiveness curve of Figure 6B in a manner described 
below. 

In a variant, some other selection may be preferred, 
e.g. from the models Fl to F3 given below, particularly 
30' -depending on the' 'type ' of nucleic' acid' that'" is"'tb be' 
quantified. 




Fl 
F2 
F3 



En = exp(-p (n-CEEP+1) )-l 

En = exp (-p (n-CEEP+1)'*) -1 

En = a - exp (-p (n-CnEP+l)°^) ~1 
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Preferably, several sets of parameters are estimated 
in step c) for several candidate changeover cycles Cg^p, 
and the minimum candidate cycle is selected for which the 
associated parameters maximize the statistical 
correlations that can be undertaken in step c) , for each 
changeover cycle Cggp. 

As mentioned above, expression (1) may also be 
written in the form: 

= fld + E,). for n e {i,2,...,n} (4) 



10 Thus, by introducing the expression (3) for 

effectiveness into formula (4), a new model is obtained 

having four parameters (Fq, Eq, p, Cggp) for the adjusted 
emitted fluorescence F„: 

fPod + E,)" f or 1 ^ n ^ C^EP 

15 = j c , ^-^^Pt-Pt"-gEep)) (5) 

iFod + E,) " "'^P*'^'-' for C,,p < n £ N 



The initial size of the population of interest, 
the effectiveness Eq of the reaction of n=0, the parameter 
P, and the changeover cycle C^j.p are evaluated 

20 repetitively for several cycle values in the vicinity of 
the changeover region CHO in order to find a statistical 
correlation maximum that is achieved for a minimum cycle 
value that is equal to the changeover cycle Cg-gp. 

In this second implementation, it is preferred to 

25 model variation in the measured and adjusted quantities 
of fluorescence as a function of cycle number on the 
basis of the models or variation in effectiveness, and 
- .subsequently -- to --carry -out - the--eorrel-at"ions —directly on" 
the measured and adjusted quantities of fluorescence. 

30 * It should be observed that by adjusting the measured 

emitted fluorescence for background noise, an artificial 
adjustment is also made on the initial fluorescence Fq. 
Thus, estimating the parameters of the effectiveness 



wo 2006/048337 



22 



PCT/EP2005/013021 



model on the basis of effectiveness 'measurements that are 
deduced from adjusted fluorescence measurements 
constitutes an additional source of error and it might be 
preferable to proceed in two stages as described below 
5 for the third implementation. 

Nevertheless, the second implementation as described 
is simpler and adapts well to PGR quantification using 
fluorescence measurements. It is based on the real 
measurements pf fluorescence F'^ which correspond to the 
10 fluorescence measurements adjusted for drift in 
background noise together with compensation e on said 
measurements. Once the background noise has been 
subtracted, we have a relationship of the following type: 

15 where 8 is a quantity that may or may not depend on cycle 
number n. It is preferably selected to be constant. 

Under such circumstances, the measured and 
"adjusted" effectiveness also written on cycle n is 

defined by: 

20 e' =11^ ^ 1 = ^"^^ ^ - 1 for all n e {l, 2,...,N ~l} - (7) 
The model of above relationship (5) thus becomes: 



F^d + s;)" f or 1 ^ n < C.^p 

^Fpd + Eo) " for Ceep < n < N 



25 Under such circumstances, the effectiveness values 

E'n approximated experimentally from the measurements 

so as to be able to set a minimum acceptable 
..effectiveness . threshold -during- >t-he- - stage of- decreasing 
effectiveness. A threshold cycle is thus determined 

30 beyond which the adjusted fluorescence measurements are 
not used for the purposes of the model (points NEG in 
Figure 4B) , . Typically, the threshold cycle corresponds 
to the first cycle in the stage of decreasing 
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effectiveness at which effectiveness drops below some 
minimum acceptable effectiveness threshold (e.g. O.IEq). 

More generally, the value of the effectiveness 
threshold preferably lies in the range 0 to 0.5, and PGR 
5 having an effectiveness value below said threshold is 
potentially biased by uncontrolled inhibition phenomena. 

In the example shown in Figure 1, the threshold 
value for effectiveness was set at 0.02 (i.e. 2% of Eq) . 
The threshold cycle . Cg corresponded to cycle n=36. 
10 Figure 7 shows the ' adjusted measurements of emitted 
fluorescence. It can be seen that there is. satisfactory 
correlation with the model (continuous line) for those 
experimental measurements (marked with an "o") up to 
cycle n=36.' Figure 8 also shows good correlation with 
15 experimental measurements for predictive effectiveness as 
obtained from Figure 7 using the model based on measured 
and adjusted fluorescence. 

The main steps of this implementation can be 
summarized as follows, with reference to Figure 9. 
20 In a start step 70, the measured values for 

quantities of fluorescence have been obtained and 
adjusted relative to background noise as a function of 
cycle number n, as shown in Figure 6A. 

In step 71, an approximation for effectiveness of 
25 the reaction in cycle n is calculated using above formula 
(2) for each of the cycles n = 1, 2, (N-1). 

In step 72, the minimum cycle Cg is determined for 
which the following two conditions are satisfied: 

the cycle Cg lies in the stage of decreasing 
30 effectiveness; and 

• the effectiveness of the threshqld cycle.. is less, 
than the threshold effectiveness value Eg (e.g. Eg = 
O.IEq): 
Ecs ^ Eg ■ 

35 It is already possible to eliminate the points NEG 

for which effectiveness is less than Eg. 



wo 2006/048337 



PCT/EP2005/013021 



24 



In step 73, a model is formed for the curve of 
adjusted emitted fluorescence which effectiveness is 
decreasing over the cycle range C^^^ = (C3 - 5) to Cg, 
using expression (8) in which it is assumed that 
5. compensation s is given by 8 = F'q: 

Thereafter, test . 74 on the value E'q estimated for 
the value E'q and the decrementation in step 75 of the 

10 value for the changeover cycle Ceep seeks to find the ' 
looked-for value of Cggp using a step size P (which may be 
equal to 1), and in repeating step 73 so long as the 
value of E'jj is less than 1. 

Thereafter, when the estimated effectiveness value 

15 exceeds the value 1 (arrow n on exiting the test 74), the 
value of the index Ceep is incremented by a step of size h 
(which may be a fraction smaller than unity) in step 76 
and in step 77 fluorescence is modeled in the same 
manner as in step 73. So long as the estimated 

20 effectiveness E'o greater than or equal to 1 in step 
78, steps 76 to 78 are repeated. When the estimated 
effectiveness takes a value of less than 1, the estimated 
parameters {F\r^\,^\,C^) are conserved in an end step 
79. 

25 In this step, a value Fq has finally been obtained 

that alone is representative of the initial population 
size Nq in the sample- of interest. It is then possible to 
use at least one standard sample having a known 
population size Nq^^ so as to determine in step 80 the 

30 initial population size Nq in the sample of interest_. 

For this purpose, a measured value of an initial 
population size F^^^ in a standard sample of known initial 
population Nost is obtained. Thereafter, the value of the 
initial population size Nq in the sample of interest is 

35 obtained by deriving a proportionality relationship 
between the measurement for the standard sample and its 
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known initial population size, and applying that 
relationship to the measurement F'q to obtain the actual 
initial population size Nq. 

In other words, in step 80 of Figure 9, it is 
5 possible .to determine the value Nq of the initial 
population size in the sample of interest by applying a 
simple proportionality relationship of the type: 
No = n (No,/Fo,) 

implying that the initial population size in the standard 

10 N^g^. and the ratio of the corrected fluorescences as 
compensated and estimated by adjusting the fluorescence 
model apply both to the sample of interest and to the 
standard sample. 

It will thus be understood that a single standard 

15 ought to be sufficient for determining the initial size 
of the population of interest in the sample of interest, 
which is an advantage provided by the invention. 

Nevertheless, in a variant, _ and where necessary, 
provision could also be made to obtain respective 

20 measured values for initial population sizes Fq^^ in a 
plurality of standard samples having known initial 
population sizes ^ost- Thereafter, a dependency 

relationship is established between the initial 
population sizes Nggt of the standard samples and the 

25 respective measured values for. their initial population 
sizes F'o5^ . Thereafter, after finding the measured value 
for the initial population size of the sample of interest 
Fq, the actual initial population size Nq of interest is 
determined by interpolation using the dependency 

30 relationship. It will be understood that this dependency 
relationship- may lalso. -typically- -be- -a- -regress-ion • of- the 
type shown in Figure 5, but having the initial 
fluorescence values Fq^j and Fq of the standards and of 
the sample of interest plotted up the ordinate (or the 

35 values of their respective logarithms) instead of 
plotting values for the changeover index C^gp. 



PCT/EP2005/013021 

26 

Once use is made of one or more standards, provision 
can be made for one or more standard samples having 
respective known initial population sizes Npst to which 
the succession of amplification reactions is applied 
5 under substantially the same conditions as for the sample 

A 

of interest. Thereafter, the measured values F'^^, for 
their initial population sizes are determined by making 
comparisons of the parameterized variations with the 
experimental measurements,, as for the sample of interest.' 

10 In other words, the same calculations are naturally 

applied concerning the measured and adjusted quantities 
of fluorescence both on the standard (s) and on the sample 
of interest. The quantity of fluorescence before the 

first cycle is estimated for the standard (s) using the 

15 same method as is used for determining F'q for the sample 
of interest, as described above. 

A third implementation, corresponding to a variant 
of the above-described second implementation consists 
overall in adjusting the model for the effectiveness E,, 

20 relative to the experimental measurements, . and in 
subsequently injecting • said adjusted effectiveness model 
into the model for the current population size N^,' or into 
the model for the measurement F^- This third 

implementation can be summari.zed as follows. 

25 The parameterized variation constructed in step b) 

is representative of yield, and in step c) , experimental 
variation of the yield is determined on the basis of' 
experimental measurements in order to compare the 
parameterized variations with the experimental variation. 

30 Thereafter, in order to obtain a parameter representative 
of the initial population size Nq the following steps are 
performed in step d) : 

dl) determining a second parameterized variation 
representative of the current population size \ in the 

35 sample of interest, making use at least of the parameter 
representing- the changeover index C^^p, and a parameter 
representative of the initial population size Nq; 
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c[2) ' applying- to said second variation^ ^ 
parameterized value for the changeover index Cg^p as 
determined in step c) / and 

d3) adjusting at least the parameter representative 
5 of the initial population size Nq by direct comparison of 
the second variation with the experimental measurements. 
Advantageously the following are performed: 

in step d2) r applying a coarse value for th& 
changeover index Cggp in the same manner as described foir 
10 detecting it with reference to above Figure 4B; and 

• in' step d3) , subsequently refining the index value 
together with adjusting the parameter representative of 
the initial population size Nq. 

Finally, it should be understood that the presently 
15 preferred second implementation of Figures 7 and 8 
differs from this third implementation by the fact that 
no attempt is made to perform correlations on 
effectiveness, but use is made merely of the mathematical 
model for effectiveness variation in order to model and 
20 refine the estimate of the correct and compensated 
fluorescence. 

Naturally, the present invention is not limited to 
the embodiments described above by way of example, and it 
extends to other variants. 

25 Thus, it will be understood that the present 

invention can also apply to relative quantification, in 
particular by PGR. In this application, as well as 
amplifying the population of interest, a reference 
population is also amplified either simultaneously in the 

30 same medium, or separately. Measurements are taken as 
follows: 

experimental measurements representative of the 
size of the population of interest; and 

experimental measurements representative of the 
•35 reference population size. 

The method ^can then continue by applying steps a) , 
b) , and c) to the reference population while step d) 
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consists merely in determining a ratio between the 
respective initial sizes of the population of interest 
and of the reference population. 

Relative quantification can be used for analyzing 
5 the expression of a gene of interest during the 
development of an organism. In order to correct in 
particular for variations in quantity and in quality 
between samples taken from the organism at different 
times, in addition to analyzing the target gene of 
10 interest, a reference gene is also analyzed that is knovn 
for having a level of expression that remains stabie 
during development. 

A final step then consists in comparing the ratios 

^Otarget 
Noref 

15 between the various samples that have been taken. 

In order to achieve the desired results, two 
strategies are possible. 

The prior art strategy is based on detecting the 
threshold cycle Ct and it normally takes place as 
20 follows. For each sample taken at different instants tO, 
tl, t2, tn, the ratio 

^Otarget 
Noref 

is determined, making use of at least one standard (i-e. 
a standard for which Notarget ^oref known) , whioh 

25 amounts to performing two successive absolute 
quantifications followed by calculating a ratio. 

Another strategy that is particularly advantageous 
in the context of the invention consists in determinirxg 
for each sample taken at different instants tO, tl, t2 , 
3*0 . • - , tn the ratio: 



\ 

^Otarget 



V Npref ^sample 
^^0 target 
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directly by using the following formula: 



N 



^ ^Otarget ^sample 
sample _ ( ^Otarget )sample_tQ 



^ ^Otarget ) ( ^Ore£ ^sample 

sample^tO 



^ ^Oref JsamplatO < ^Oref )sample^tO 



In this second implementation, which in the end 
5 makes use only, of the parameter Fq, in combination with 
the technique of the invention, no standard sample is 
needed, which is particularly advantageous. 

Reference is now made to Figure 10 which shows an 
installation for implementing the method of the 

10 invention. It comprises a support SUPP in this case 
comprising a well containing the sample of interest ECH 
and a well containing a standard sample referenced St, 
for example. The support SUPP * is enclosed in an 
enclosure ENC, e.g. fitted with heater means (not shown) 

15 for applying a PGR reaction to the standard and to the 
sample of interest. 

In the example described, provision is preferably 
made to take measurements of the quantities of 
fluorescence emitted on each cycle, both by the standard 

20 St and by the sample of interest ECH. To this end, a 
selected reagent is inserted into the wells and the 
samples are illuminated by a lamp (e.g. a halogen- 
tungsten lamp) in order to measure the respective 
quantities of fluorescence coming from the sample of 

25 interest and from the standard sample on each PGR cycle 
that is applied. In addition, an apparatus for detecting 
fluorescence comprises, for example, an objective lens 11 
for ''coirec:ting'*'the * light"" coming"'* from " the fluorescence, 
and photon counting means 10, e.g. a charge-coupled 

30 device (GGD) camera, and/or photomultipliers, in order to 
measure the fluorescence emitted on each PGR cycle from 
the sample of interest and from the standard. Thus, the 
fluorescence emitted by each well is advantageously 
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focused by the lens" 11 and then is preferably detected by 
a CCD camera 10 connected to an acquisition card 21, e.g. 
of the Personal Computer Memory Card Internatiional 
Association (PCMCIA) type provided in a central unit 20 
5 of a computer. 

The computer is then connected to the above- 
mentioned measuring means 10 to receive therefrom signals 
that are representative of the measured quantities of 
fluorescence detected on each PCR cycle, and to pr-ocess 
10 these signals in order to determine an initial size for 
the population of interest prior to the first cycle, by 
implementing the method of the invention. 

Typically, the processor unit- comprises the 
following: 

15 . an acquisition card 21 connected to the 

measurement means 10; 

■ working memory 25 (e.g. of the random access 
memory (RAM) type) for temporary storage and processing 
of the above-mentioned signals; 

20 • permanent memory 24 for storing the computer 

program product in the meaning of the invention and for 
storing the data that has been processed and that is 
ready for use, e.g. in subsequent diagnosis; 

• where appropriate, a reader 22 of a memory medium 
25 such as a CD-ROM, for example, which may initially carry 

the computer program product; 

optionally a communications interface 2 6 for 
communicating with a local or remote site (connection 
28), e.g. for transmitting the processed data so as to 
30 enable a diagnosis to be made remotely concerning a 
patient; 

• a /graphics "i^ 27 typically connected to a 
display screen 30; and 

a processor 23 for ■ managing the- interactiions 
35 between these various items of equipment. 
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The computer may also have input members such as a 
keyboard 41 and/or a mouse 42 connected to the central 
unit 20 • 

Nevertheless, it should be understood that in tine 
5 meaning of the invention the installation comprises 
overall: 

• a sample support SUPP, at least for the sample of 
interest; 

* a first apparatus ENC for applying said succession 
10 of amplification reactions at least to the population of 

interest in the sample of interest; 

a second apparatus 10 for taking measurement's 
representative of the current size of the population of 
interest; and 

15 . . computer means 20 suitable for receiving 

measurement signals from the second apparatus 10 and for 
implementing all or some of the steps of the method of 
the invention. 

For this purpose, a computer program product can toe 

20 used for controlling the computer means. The program may 
be stored in a memory of the processor unit 20 or on a 
removable memory medium (CD-ROM etc.) and suitable for 
co-operating with the reader of the processor unit. Tb.e 
computer program in the meaning of the invention th&n 

25 contains instructions for implementing all or some of th-e 
steps of the method of the invention. For example, th-e 
algorithm of the program may be represented by a flo-w 
chart equivalent to the diagram of Figure 9. 



